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Abstract 

On the basis of a semi-phenomenological model, it is argued that the high energy anomaly 
observed in recent photoemission experiments on cuprates is caused by interaction with an over- 
damped bosonic mode in the mid-infrared region of the spectrum. Analysis of optical conductivity 
allows to connect this excitation to the incoherent charge response reported for the majority of 
high Tc materials and some other perovskites. We show that its large damping is an essential 
feature responsible for the "waterfall" dispersion and linewidth of the spectral weight. 

PACS numbers: 74.72.-h,79.60.-i, 74.25.Jb 
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The "high energy anomaly" (HEA) in the angle resolved photoemission spectra of hole 
and electron doped cuprates has been reported recently in a number of works [ll-lsl. The 
anomaly is characterized by the presence of high energy kink in the dispersion derived from 
angle resolved photoemission spectroscopy (ARPES) at about 0.2 — 0.4 eV below the Fermi 
surface followed by a puzzling "waterfall" structure at larger energies. The latter is highly 
incoherent and follows an almost vertical dispersion around the F— point extending till 0.9 
eV. When measured along the main symmetry directions of the Brillouin zone such broad 
linewidth remains roughly constant along the waterfall |l, The momentum distribution 
of the spectral weight at a fixed energy above the kink resembles a diamond shape (e.g. for 
E = 0.4 in BISCO 2212 with the corner at K ~ (vr/S, 0) and side center at ~ (7r/4, 7r/4)). 
At even higher energies the spectrum recovers the parabolic dispersion reminiscent of the 
"bare" band. Several mechanisms responsible for the observed behavior have been proposed, 
such as disintegration of the quasiparticles into a spinon and holon [l| , polaronic effects Ml , 
spin fluctuations whose spectrum can reach the relevant energy scale 13, l3, H^, coherence- 

nnn 

incoherence crossover [5|, la, 19|. However a controversy also exists in experimental data, e.g. 
regarding the role of matrix element effects in ARPES , luj ) . The fixed energy slices 

of the incoherent tails reveal the existence of a "grid-like" structure when the spectral weight 
maxima are aligned to the (0, 1) , (1, 0) directions and suppressed in the diagonal directions. 
It was suggested that selfenergy effects may therefore not be a cause of the anomaly which 
instead should be related to the presence of one dimensional structure. On the other side 
some experiments on magneto-optics spanning a broad frequency window (4000 cm~^) are 
apparently inconsistent with magnetic excitation scenario 121. Il3|. 

In Ref. jl^ the optical conductivity spectra have been analysed by a method allowing to 
extract the electron selfenergy within apparently general assumptions about the spectrum. 
However, no indication of the HEA has been found, while the well known low energy kink 
(LEK) and the above mentioned feature at 0.9 eV are reproduced in good agreement with 
ARPES. On the other side, the analysis based on fermion-boson concept in Ref. ISj before 
the discovery of HEA has revealed that optical data do indicate the existence of a high 
energy cutoff of the boson spectrum at around 300 meV. Also in several other earlier papers 
treating electron-boson scenario it has been demonstrated that taking into account the 
: inite linewidth of the boson is important for the faithful description of the ARPES spectra 



id llTI . Il8| . In particular, the LEK is well reproduced if the linewidth of the boson (phonon) 
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is about 5 — 10% of its energy. In this paper we show that the quahtatively new features of the 
HEA can be accounted for by assuming the "overdamped" regime for the "boson", i.e., that 
the hnewidth of this excitation is significantly larger than the characteristic energy. This 
assumption is motivated by the ubiquity of such incoherent excitations in the mid-infrared 
(MIR) frequency range of the optical spectra in cuprates as revealed, e.g., by Drude-Lorentz 
analysis. As explained below, although the bare frequency of such excitation is larger than 
the high energy kink, the scale that actually characterizes its effect on electron excitations 
falls in the relevant energy range. We consider the selfenergy effects within a simple model of 
electrons interacting with a bosonic mode in the form of damped oscillator. The parameters 
of the oscillator are chosen to match the known experimental data, in particular those from 
ARPES. As for instance, its coupling as extracted from the slope of single particle spectrum 
in the region preceding the HEA (e.g. for BISCO 2212 and LBCO is Ac = 1 — 1.5), momentum 
and energy location of the observed structures etc. [s], lol. Calculations show that the 
constraints set by the model and experimental data are rather restrictive so that the choice 
of parameters can not be reduced to simple fitting. For instance, it is not possible to obtain 
a vertical dispersion by choosing a small damping parameter and varying the others. We 
show that, indeed, the key component of the oscillator that accounts for the main anomalies 
like the "waterfall" and the large width of the ARPES tails, is the damping parameter. 
We further find quantitative similarities between the parameters of the oscillator and data 
available from optical experiments on different high Tc materials, namely, the incoherent 
excitation spectrum in the MIR region. It is known that this mid-gap incoherent part of 
;he spectral weight emerges more quickly than the "coherent" Drude part upon doping 
2o| . Unfortunately, optical data can be related only to the long wave limit {q = 0) of the 
"overdamped boson" , as there are no experimental probes available to study its dispersion 
characteristics. Therefore for our model calculations we have neglected the dispersion of this 
MIR oscillator. Nevertheless, we discuss some features reported in experiments that can be 
accounted for by considering an anisotropic dispersion. 

In the Table I we have summarized the experimental data for the Lorentz oscillator 
component in the relevant MIR region of energies extracted from optical conductivity for 
the Drude-Lorentz representation of the dielectric function 




(1) 
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The Drude component with damping 71 and weight determined by the plasmon frequency 
Up^i (~ 1.5 — 2 eV) is centered at u; = 0, the frequency and damping of the discussed 
excitation are further denoted by Uq and 7, the "effective" frequency Q is explained below. 
The values of the oscillator strength corresponding to ujq can be found in the listed papers, 
but are generally comparable to the weight of the Drude part of the dielectric function at 
optimal doping. It should be noted that data in the table are both doping and temperature 
dependent and that respective multi-component analysis in terms Drude-Lorentz oscillators 
depends on the number of components chosen. In particular, physical arguments suggest 
that an oscillator in the phononic frequency range of the MIR should be distinguished 
from the oscillator singled out in the table. Such a fitting procedure can be formalized. 



for instance, in the Kramers-Kronig constrained variational analysis, e.g. 
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28|. The 



obvious deficiency of such analysis is the absence of frequency dependence in the mass 



and relaxation rate of the Drude component (generalized Drude model) [29|. However for 
the energies of interest in the present context it is not expected to dramatically alter the 
estimations. Indeed, the generalized Drude model is appropriate for the lower energy part of 
the spectrum characterized by a single scattering mechanism, but is insufficient to account 
for the higher energy scale related to the HEA. For the estimates discussed in the present 
paper we assume that these scales can be treated independently in calculating the electronic 
selfenergy. As already mentioned, the long wavelength values of characteristic frequencies 
in the Table should be viewed as an estimate for the actual frequency of the overdamped 
charge excitation for our model since its value at finite momentum is not known. 

We introduce the retarded bosonic propagator B (g, u) in the form of damped oscillator 
jsol with spectral density ps {q,^^)'- 

The choice of ([2]) in our model is motivated by its compliance with the generic form of the 
response functions frequently used in analyzing experimental data that, in turn, corresponds 
to physical requirements, such as Kramers-Kronig relations. Coupling to electron excitations 
is considered in the second order of perturbation theory, i.e. non-selfconsistently. The 



scale of the boson frequency is assumed to still a.. 



reasonable at least for evaluation purposes (see |31i for a recent discussion). Although the 



low the adiabatic approximation to be 



vertex corrections and selfconsistent treatment are beyond the scope of the paper, their 



analysis certainly deserves a special consideration. Using the spectral representation for 
the Matsubara Green function, after carrying out the frequency summation and analytic 
continuation to real frequencies one obtains the second order expression for the retarded 



electron selfenergy, e.g. 



18|, 



32| 



ub (y) + riF (x) 



POO POO 

T.R{k,uj) = — S^g'^ {k,q) j / pp {k - q, x) ps iq,y) ' '""^ dxdy, (3) 

where g (q) is the coupling constant and ub^f are Boson, Fermion distribution functions. 
The linewidth of the electronic spectral function pp {k — q,x) is neglected, pp {k — q,x) = 
6 {x — Sk-q) 1 and a tight binding parametrization of the band dispersion is used as explained 
above, that allows to focus on the scattering mechanism due to the "boson". If also the 
linewidth of (^^) is neglected, the expression ([3]) transforms into the one that has been 
often employed in describing the LEK ( see, e.g.[33|) and results in the usual logarithmic 
singularity of ReS {uj) at the boson frequency. To study the effect of the broad linewidth 
on electron spectrum we neglect the dispersion in spectral function ([2]) and in the coupling 
g (g) . Then the selfenergy becomes a local function 



oo 



(.) ^ X.,.l I dy^y^, [ dEp, iE) ( ,--/^^+:-(fl ) , (4) 



where po (E) is the electronic density of states, Ac is the dimensionless coupling constant. 
The coupling constant A determined from the slope of Re S {u) is slightly different from Ac. 
The band width of Ek is taken as the energy unit in numerical calculations. 

The spectral density of the damped oscillator ([2]) is peaked at the frequency Q = 

1/2 



(12 - 6i?2 + 6V16 - 4i?2 + i?4j X uo/Q, where R = 'j/luo. It is clear that the effective 
frequency fl is significantly lower than uo due to a large damping typical for the MIR os- 
cillator. As seen from the FigUl it is this effective frequency that sets the location of the 
high energy kink. In the figure we plot the ReS (uj) and ImS {u) together with the "boson" 
spectral function ps (uj) to make the correlation between the peaks clearly visible. The slope 
of the ReS {uj) corresponds to the coupling constant 1.3 and the ratio of t' /t = 0.1 was taken 
to model the LSCO band Ek- One can also see that the {—ImH {uj)) follows a tu^ at low 
energies and its inflection point corresponds to Vt. After reaching its maximum at higher 
energy {—ImH {uj)) follows an almost fiat region until it starts decreasing closer to the band 
bottom. Such behavior has been reported in l|, |5| for the linewidth of the ARPES spectral 



function. In our calculations we clearly see the feature mentioned in 



IJ] and related to 



the finite bandwidth when the renormalized dispersion curve crosses with bare one. The 
crossing occurs at certain nonzero frequency within the waterfall region where ReS {uj) = 0. 

The calculated dispersion of the spectral function peak along the two main symmetry 
directions of the Brillouin zone are presented in Figl2]for doping 5 = 0.15 and band structure 
of BISCO 2212. The dashed curve shows the evolution of the curve at larger doping, 6 = 0.25, 
for the diagonal direction. The HEA is close to ko ^ tt/A at i?i ^ i7 = 0.15 (~ 0.4 eV), the 
"waterfall" is superceded by parabolic dispersion at E2 — 0.33 0.85 eV). The momentum 
of HEA ko decreases with doping. These results qualitatively agree with experimental data 
[1, sl and appear to suggest that the band structure is of primary importance for the location 
and doping dependence of the HEA, as the doping dependence of /cq follows the shrinking 
of the Fermi surface. However, it should be taken into account that also the parameters of 
the MIR oscillator depend on doping. In particular, the quantitatively weaker dependence 
of ko observed in experiments could be related to the decrease of uq with doping. We are 
not aware of the doping dependence of the linewidth of the spectral function in the waterfall 
region, but our calculation (not shown) indicates its decrease with doping. The temperature 
dependence of S (u) is very weak in the relevant region of energies. We find a broader 
linewidth in the (0, 1) and (1,0) directions as compared to the diagonal in agreement with 
the measurements presented in ^]. Figj3] shows the calculated three dimensional dispersion 
of the spectral density peak within the Brillouin zone obtained byneglecting the dispersion 
of the damped oscillator that agrees with the results reported in [l| . However, as noted, the 



ARPES measurements presented in some other papers [3|, |ll| show an anisotropic structure 
of the HEA. These results can be described within the present model by introducing the 
anisotropy into the parameters of the oscillator. For instance, the difference in energy of the 
kink for momentum cuts along different directions of the Brilouin zone found in d\ can be 
explained by a larger value of Q along the diagonal direction. In a different context, it has 
been noted that the MIR peak in optics is affected by lattice oscillations 3^ . It has also been 
suggested 35| that the anomalous softening of the bond-stretching in plane phonon mode 
is related to an overdamped charge excitation in the same MIR part of the spectrum and 
that the behavior of this mode could serve as an energy-momentum resolved probe for the 
charge excitation. This raises an intriguing possibility that the same overdamped excitation 
is responsible for the phonon anomaly and the HEA. One of its characteristics implied from 
the phonon analysis is the anisotropy of the charge response at finite momentum in the 
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energy range that matches the grid-hke pattern mentioned in ll|. Namely that anomalous 
phonon softening occurs in the (0, 1) , (1, 0) directions and a much smaller effect is observed 
in the (1, 1) direction. This anisotropy is related to a larger charge response located at lower 
energies for the orthogonal direction as compared to more dispersive and weaker response 
along the diagonal. The strength of the anisotropy depends on the material and doping. 
Interestingly, in the papers 0, sl the energy of the kink (in absolute value) in LBCO, 
BISCO 2212 was found to decrease significantly away from the diagonal direction of the 
Brillouin zone. This behavior is also in agreement with the MIR charge "boson" scenario 
proposed here. The broad incoherent spectrum above the high energy kink emerges due 
to the composite nature of electronic excitation (in itself a source of finite linewidth) with 
emission of an "overdamped boson". This explains the qualitative difference between the 
high and low energy kinks. In the latter case the linewidth of the boson is not more than 10% 
of its energy scale, while in the case of MIR it is an order of magnitude larger. At the same 
time, for a comparable coupling constant, such linewidth accounts for the verticality of the 
dispersion and linewidth of the waterfall. As we have seen, the ImS remains almost fiat in 
this region. At energies around 0.9 eV the ImE starts to decrease and the dispersion recovers 
the bare band parabolic shape, albeit shifted to higher energies. It should be mentioned that 
the optical conductivity data have served only for estimations and, as discussed above, one 
needs to have a probe for momentum resolved charge excitations in the 0.4 — 0.5 eV range 
to characterize the "overdamped boson". 

In |l| several other types of materials have been scrutinised for HEA, but no evidence 
has been found. It follows from the proposed scenario that HEA should exist in the ARPES 
of other perovskite materials where the incoherent MIR excitation has been observed (and 
vice versa). For instance, in the paper {s^ the ARPES analysis is focused on the LEK at 60 
meV in bilayer manganite La2-2xSri^2xMn20i. However, also the data for higher energies 
are presented, where the existence of the other kink at around 0.45 eV is clearly visible. It 
would also be interesting to carry out a similar search for the HEA in ARPES experiments 
on a non-magnetic perovskite such as Bai^xK^BiOs with a relatively high Tc that is also 
known to have an incoherent charge excitation in the MIR about 0.4 — 0.6 eV 



33, 



38 



39| 
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as well as anomalous softening of the bond-stretching phonon mode 



40| 
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FIG. 1: Real and imaginary parts of the selfenergy S (w) (jl]) resulting from interaction with 
the overdamped oscillator with spectral density pb{uj). Parameters are lvq = 0.18,7 = 0.3, A = 
1.3, T = 0.001, 5 = 0.15. The flattening of ImS (uj) corresponds to the incoherent "waterfall" region 
in ARPES with almost constant linewidth. 
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FIG. 2: Dispersion of the spectral density maximum along the diagonal direction, K = Kx = Ky, 
(continuous line 6 = 0.15, dashed line 6 = 0.25) and orthogonal direction, K = Kx{Ky), Ky{Kx) = 
0, (dotted line 5 = 0.15) of the Brillouin zone. Parameters are: loq = 0.18,7 = 0.21, A = 1.3, T = 
0.001. 




FIG. 3: Three dimensional plot of the surface defining the position of the peak in the spectral 
density as a function of energy and in-plane momentum. Parameters correspond to FiglJJ 
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TABLE I: Parameters of the MIR oscillator in ([T]) determined from the Drude-Lorentz analysis of 
the optical conductivity. 



material 
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VI 


Quiiada.[21] Bi2Sr2CaCu20s,a. axis 


0.545 


0.95 


0.358 


Quiiada.[211 Bi2Sr2CaCu208,h axis 


0.51 


1 


0.3 


Wang, [22] Lai,gCai,iCu20Q+s 


0.59 


0.74 


0.47 


Wang,[22] Lai,85Sro,i5CaCu206+s 


0.57 


0.67 


0.47 


Uchida, [23] Lai^gSro^iCuO^-^s 


0.46 


0.8 


0.3 


Kircher,[24] YBa2Cu408 


0.55 


1.12 


0.31 


Heumen,[25] HgBa2Cu04^s 


0.48 


0.9 


0.3 


Carbone, [26] Bi2Sr2Ca2Cu-iOiQ 


0.63 


1 


0.44 


Orenstein, [27] Y Ba2Cu^OQ,g 


0.65 


1.3 


0.38 
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